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The production of J/ip mesons in continuum e + e~ annihilations has been studied with the BABAR 
detector at energies near the T(4S) resonance. The mesons are distinguished from J/vj production 
in B decays through their center-of-mass momentum and energy. We measure the cross section 
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e + e~ -> J/ipX to be 2.52 ± 0.21 ± 0.21 pb. We set a 90% CL upper limit on the branching fraction 
for direct T(45) -> J/tpX decays at 4.7 x 10" 4 . 

PACS numbers: 13.65.+i, 13.25.Gv, 12.38.Qk, 14.40.Gx 



The development of non-relativistic QCD (NRQCD) 
represents a significant advance in the theory of the pro- 
duction of heavy quarkonium (qq) states In partic- 
ular, it provides an explanation |2| for the cross section 
for ip(2S) production observed by CDF which is a 
factor of 30 larger than expected from previous models. 
The enhancement is attributed to the production of a cc 
pair in a color octet state, which then evolves into the 
charmonium (cc) meson along with other light hadrons. 
A similar contribution is expected in NRQCD for J/tp 
production in e + e~ annihilation 0,0, but is absent in 
the color singlet model [6|. 

Significant continuum J/ip production — as distinct 
from production in B decay at the T(4S) resonance — has 
not been observed previously in e + e~ annihilation below 
the Z resonance. It therefore represents a good test of 
NRQCD. In particular, matrix elements extracted from 
different J/ip production processes should be consistent 
0- In addition, momentum, polarization and particu- 
larly the angular distributions of the J/ip distinguish 
between theoretical approaches ||. Despite NRQCD's 
successes, it is not clear that it correctly explains [9j the 
CDF measurements of J/ip polarization I1C|. or measure- 
ments of J/ip photoproduction at HERA[lll ll2T ] . 

The study reported here uses 20.7fb~ 1 of data col- 
lected at the T(45) resonance (10.58 GeV) and 2.59 fb" 1 
collected at 10.54 GeV, below the threshold for BB cre- 
ation. The luminosity-weighted center-of-mass (CM) en- 
ergy is 10.57 GeV. 

The data were collected with the BABAR detector 01 
located at the PEP-II collider at the Stanford Linear Ac- 
celerator Center. PEP-II collides 9 GeV electrons with 
3.1 GeV positrons to create a center of mass moving along 
the z axis with a Lorentz boost of 0^ — 0.56. 

The momenta and trajectories of charged particles are 
reconstructed with two detector systems located in a 1.5- 
T solenoidal magnetic field: a five-layer, double-sided sil- 
icon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH). The fiducial volume covers the polar angular re- 
gion 0.41 < 9 < 2.54 rad, 86% of the solid angle in the 
CM frame. 

The energies of electrons and photons are measured in 
a CsI(Tl) electromagnetic calorimeter (EMC) in the fidu- 
cial volume 0.41 < 9 < 2.41 rad, 84% of the solid angle 
in the CM frame. The instrumented flux return (IFR) is 
used to detect muons. The DIRC, a unique Cherenkov 
radiation detection device, distinguishes charged parti- 
cles of different masses. 

J/ip mesons are reconstructed via decays to electron or 
muon pairs. The leptons must form high-quality tracks 
with 0.41 < 9 < 2.41 rad: they must have p t > 0.1 GeV/c 



and momentum below 10 GeV/c, have at least 12 hits in 
the DCH, and approach within 10 cm of the beam spot 
in z and within 1.5 cm of the beam line. The beam spot 
rms size is approximately 0.9 cm in z, 120 /im horizontally 
and 5.6 /im vertically. 

One electron candidate must have an energy deposit 
in the EMC of at least 75% of its momentum. The other 
must have between 89% and 120%, and must also have an 
energy deposition in the DCH and a signal in the DIRC 
consistent with expectations for an electron. Both must 
satisfy criteria on the shape of the EMC deposit. If pos- 
sible, photons radiated by electrons traversing material 
prior to the DCH are combined with the track. 

Muon candidates must deposit less than 0.5 GeV in the 
EMC (2.3 times the minimum-ionizing peak), penetrate 
at least two interaction lengths A of material, and have a 
pattern of hits consistent with the trajectory of a muon. 
We require that the material traversed by one candidate 
be within 1 A of that expected for a muon; for the other 
candidate, this is relaxed to 2 A. 

The mass of the J/ip candidate is calculated after con- 
straining the two lepton candidates to a common origin. 

To reject interactions with residual gas in the beam 
pipe or with the beam pipe wall, we construct an event 
vertex using all tracks in the fiducial volume and require 
it to be located within 6 cm of the beam spot in z and 
within 0.5 cm of the beam line. To suppress a substan- 
tial background from radiative Bhabha (e + e~7) events 
in which the photon converts to an e + e _ pair, five tracks 
are required in events with a J/ip — + e + e~ candidate. 

At this point, the data includes J/ip mesons both from 
our signal — continuum-produced Jji\j mesons and J/ip 
mesons from the decay of continuum- produced ip(2S) 
and Xcj mesons — and from other known sources. We ap- 
ply additional selection criteria to suppress these other 
sources based on their kinematic properties. 

The most copious background, B — > J/0 A, is elimi- 
nated by requiring the J/ip momentum in the CM frame 
(p*) to be greater than 2GeV/c, above the kinematic 
limit for B decays. This requirement is dropped for data 
recorded below the T(AS) resonance. 

Other background sources include initial-state radia- 
tion (ISR) production of J/ip mesons, e + e~ — > "fJ/ip, or 
of the ip(2S), with ip(2S) -> J/ipX. ISR production of 
lower-mass T resonances is negligible. Two photon pro- 
duction of the %c2 can produce J/ip mesons via \ci — * 
7 J/ip . Because the out-going electron and positron are 
rarely reconstructed, this process, like the ISR J/ip pro- 
duction, contains only two tracks. We therefore require 
three high-quality tracks with 0.41 < 9 < 2.54 rad. 

The remaining background is primarily ISR ip(2S) de- 
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cays to J/ip 7r + 7r - , plus some ISR J/tp events in which 
the ISR photon converts. To suppress these, we require 
the visible energy E to be greater than 5 GeV, and the 
ratio of the second to the zeroth Fox- Wolfram moment 
(l4|. i?2, to be less than 0.5. Both are calculated from 
tracks and neutral clusters in the fiducial volume. Fig- 
ure ^ which displays the visible energy and i?2 distribu- 
tions for our signal and for simulated ISR background, 
motivates these criteria. 

The ISR distributions in Fig. ^ are obtained from a 
full detector simulation. All selection criteria are ap- 
plied, other than the one on the quantity being plotted. 
ISR kinematics ensures E < 5 GeV when the photon is 
outside the fiducial volume unless it interacts in material 
and deposits additional energy in the detector. The rate 
of such interactions is not accurately simulated and so is 
obtained by a comparison to data for E < 5 GeV. Ap- 
proximately 3.5% of the J/tp meson events that satisfy all 
criteria are from this background; an additional ~ 1.6% 
are ISR events with the photon in the fiducial volume. 
Systematic errors on the remaining backgrounds are es- 
timated from a comparison between simulation and data 
for E < 5 GeV and for events in which the ISR photon is 
reconstructed. 

J/ip production as a function of E is obtained in data 
by fitting the dilepton mass distribution in 1-GeV wide 
energy intervals after applying all other selection criteria. 
The fit uses a polynomial function for the background 
distribution. The J/ip mass function is obtained from 
a complete simulation of B — > J/tpX events, convolved 
with a Gaussian distribution to match the resolution of 
12 MeV/c 2 observed in data in a sample of approximately 
14,000 B — > J/tpX events. The signal distribution in E 
is obtained by subtracting the ISR backgrounds from the 
data distribution. 

A similar process is used for i?2- Figures ^c) and (d) 
show there is little signal above i?2 of 0.5. In this respect, 
the continuum J/ip events are more similar to BB events, 
in which the energy is distributed spherically, than cc 
events, which tend to be jet-like. 

The mass distributions of the selected J/ip candidates 
show clear signals for both e + e~ and final states, 

both on and below resonance (Fig. EJ) . 

To determine the production cross section, we perform 
mass fits in 15 p* - cos 9* bins, where 9* is the polar 
angle of the candidate in the CM frame. This allows 
us to correct for the variation of efficiency with p* and 
cos 9*. The cross section is given by: 



-^J/ipX 



(JV, 



B 



13) 



(1) 



where the sum is over three p* (i) and five cos 9* (J) 
bins. Nij is the number of J/ip mesons in the bin, where 
electrons and muons are analyzed separately, but off and 
on-resonance data are combined. The sum of the yields 
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FIG. 1: J/ip signal events observed as a function of visible 
energy E in the (a) e + e~ and (b) fj, + fi~ final states; R2 
distribution for (c) e + e~ and (d) fi + fi~ . The histogram is 
the predicted ISR background that has been subtracted from 
data; the filled histogram is the ISR ip(2S) component only. A 
requirement of > 5 tracks is applied to the e + e~ sample; ap- 
plying it to the /i + fi~ sample produces the dashed histogram. 
Event preselection requires E > 4 GeV and R2 < 0.95. 
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FIG. 2: Mass distribution of J/ip candidates reconstructed in 
data recorded below the T(4S) resonance in the (a) e + e~ and 
(b) fj, + final states. Mass distributions for p* > 2 GeV/c in 
data at the T(4S) resonance in (c) e + e~ and (d) fi , n~ final 
states. The number of J/ip mesons extracted by a fit to the 
distribution is shown on each graph. 
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from the 15 fits agrees to within 1% with the yields in 
Fig. [21 Bij is the ISR background, B.i M.^e +e- is the 



J/tp — > e + e _ or branching fraction 15], and Ci is 

the integrated luminosity — sum of on plus off-resonance 
for p* > 2 GeV/ c, off-resonance only for p* < 2 GeV/ c. 

The reconstruction efficiency (acceptance, track 
quality and lepton identification) is calculated in each bin 
with simulated unpolarized J/ip mesons uniformly dis- 
tributed in p* and cosO* . The efficiency decreases with 
increasing p* and cos 8* due to acceptance. The average 



& is 0.63 for J^j 



and 0.48 for J/tj; 



A t+ A t ~ 



where the difference is due to lepton identification. 

Particle identification efficiency is verified in data by 
comparing the number of J/tp mesons in B decays in 
which one or both leptons satisfy the requirements. The 
efficiency of the track-quality selection is studied by com- 
paring tracks found in the SVT and DCH. 

The components of e E , the event selection efficiency, 
are determined as follows. We estimate the efficiency of 
the requirements on the number of high-quality tracks, 
primary vertex location, and total energy to be the aver- 
age of simulated cc and BB events, and the uncertainty 
to be one-half the difference. We use BB events for i?2- 

The efficiency of the five track requirement applied to 
e + e _ candidates is 0.67, obtained by comparing the net 
J/ip yield, in both e + e~ and \i + [i~ final states, in events 
passing and failing the requirement. Overall, e E = 0.59 
for e + e _ and 0.89 for 

The calculations of the J/tpX cross section from the 
e + e~ and /i + ^t~ final states are consistent: the ratio 
a( f i+fi-)/a(e+e-) is 0.93 ±0.11 for p* > 2 GeV/c. The 
two values are combined, accounting for common system- 
atic errors, to obtain 



Te+e-^X = 2.52 ± 0.21 ± 0.21 pb, 



(2) 



where the first error is statistical and the second sys- 
tematic. With existing values for matrix elements, color 
singlet cross section estimates range from 0.45 to 0.81 pb 
while NRQCD cross sections, including a color 
octet component, range from 1.1 to 1.6 pb 0,O|. 

The dominant component of the 8.3% systematic error 
is a 7.2% uncertainty on e E common to both the e + e~ 
and cases and a 4.9% uncertainty due to the five 

track requirement. Other contributions include 2.4% due 
to track quality cuts; 1.5% from the luminosity; 1.8% 
(electrons) or 1.4% (muons) from particle identification; 
and 1.2% from the ISR background. 

The statistical error is dominated by the uncertainty 
on the contribution below p* of 2 GeV/c. Restricting 
the measurement to p* > 2 GeV/c gives <r e +e-^j/ipx — 
1.87 ± 0.10 ±0.15pb. 

In determining the cross sections, we assume that 
there are no J/ip mesons from direct T(AS) decays. We 
quantify this statement using the p* > 2 GeV/c com- 
ponent. We scale the off-resonance event yield to the 
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FIG. 3: Center of mass momentum distribution of J/ip mesons 
produced in continuum e + e~ annihilation. 



on-resonance luminosity and subtract it from the on- 
resonance yield. The excess, attributable to T(4S) de- 
cays, is consistent with zero: — 120±179 e + e~ events and 
176 ± 138 /Lt+M", in a sample of (22.7 ± 0.4) x 10 6 Y(4S) 
decays. Using the average reconstruction efficiency for 
p* > 2GeV/c (0.62 for e+e" and 0.45 for fi+fjT), we ob- 
tain B T{4S) ^ J/4 ,x = (1-5 ± 2.2 ± 0.1) x 10" 4 . A Bayesian 
90% confidence level upper limit with a uniform prior 
above zero is: 



Br ( 4S)^x < 4.7 x 10~ 4 (90% CL), 



(3) 



for J/ip with p* > 2 GeV/c. This result disagrees with 
a previous publication 16]. In NRQCD, the expected 
partial width is similar to that for the T(IS) ja, Il7]] . 
implying a branching fraction of a few x 10~ 6 . Note that 
a true branching fraction of 10 -4 would correspond to an 
effective cross section of 0.10 pb. 

Production and decay properties of the J/ip have also 
been studied. The p* distribution is obtained by divid- 
ing the sample into 500MeV/c wide intervals, fitting the 
resulting mass distribution, subtracting predicted ISR 
backgrounds, correcting for the reconstruction efficiency, 
and normalizing for different luminosities (Fig. [3J - 

The distribution of the signal in cos#* has been ex- 
tracted and fit with 1 + A ■ cos 2 6*. Both NRQCD 
and color singlet calculations predict a flat distribution 
(A w 0) at low p*. At high momentum, NRQCD pre- 
dicts 0.6 < A < 1.0 while the color singlet model pre- 
dicts A w —0.8 8]. We measure the distribution sep- 
arately for low and high momentum mesons, selecting 
p* = 3.5GeV/c as the boundary. We proceed as for the 
p* distribution, with mass fits performed in cos 9* inter- 
vals of width 0.4. The distributions are then normalized 
to unit area (Fig. EH). We find A = 0.05 ± 0.22 for 
p* < 3.5GeV/c and A = 1.5 ± 0.6 for p* > 3.5GeV/c, 
clearly favoring NRQCD. 

Finally, we obtain the helicity angle Or distribution 
for the two p* ranges by fitting mass distributions in in- 
tervals of width 0.4 in cosOh (Fig. 0Jd). The helicity is 
the angle, measured in the rest frame of the J/ip, be- 
tween the positively charged lepton daughter and the di- 
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FIG. 4: (a) Production angle (cos#*) distribution for J/ip 
mesons produced in continuum e + e~ annihilation; (b) helicity 
(cos#h) distribution. Solid curve is the fit to p* < 3.5 GeV/c; 
dashed curve is for p* > 3.5 GeV/c. 



rection of the J/tp measured in the CM frame. Fitting the 
function 3 (l + a ■ cos 2 Oh) /2 (a + 3), we obtain a J/tp 
polarization a = -0.46 ± 0.21 for p* < 3.5GeV/c and 
a = -0.80 ± 0.09 for p* > 3.5 GeV/c. a = indicates an 
unpolarized distribution, a — 1 transversely polarized, 
and a = — 1 longitudinally polarized. 

In summary we measure the cross section 
(Je+e-^j/ipx = 2.52 ± 0.21 ± 0.21 pb. Restricting 
to p* > 2GeV/c, we find 1.87 ± 0.10 ± 0.15 pb. The 
total cross section and the angular distribution favor 
the NRQCD calculation over the color singlet model. 
We set a 90% CL upper limit on the branching fraction 
T(45) -> J/4>X of 4.7 x 10~ 4 . 
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